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We study the effects of weak off-diagonal disorder on Peierls systems 



with a doubly degenerate ground state. We show that for these systems 
disorder in the electron hopping amplitudes induces a finite density of 
solitons in the minimal-energy lattice configuration of a single chain. 

. . i 

<^ \ These disorder- induced dimerization kinks are neutral and have spin ±. 

Using a continuum model for the Peierls chain and treating the lattice 

T 1 

Q\ \ classically, we analytically calculate the average free energy and density 

of kinks. We compare these results to numerical calculations for a dis- 
crete model and discuss the implications of the kinks for the optical and 
q magnetic properties of the conjugated polymer trans-p olyacetylene. 



PACS numbers: 71.20.Rv, 63.20.Kr, 71.23.-k 



I. INTRODUCTION 



Recently we considered the effects of weak disorder in the electron hopping amplitudes 
on the lattice configuration of Peierls systems with a doubly degenerate ground state, 
using the conjugated polymer trcms-polyacetyiene as an exampie.0 In the absence of 
disorder the ground state of a irans-polyacetylene chain is uniformly dimerized due to 
the interaction of the lattice with the half-filled band of 7r-electrons propagating along 
the chainJl The ground state energy is independent of the sign of the dimerization, i.e., 
it is equal for the two carbon-carbon bond alternations . . .-long-short-long-short-. . . and 
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. . .-short-long-short-long-. . .. Disorder in the electron hopping amplitudes, originating, 
e.g., from random twists of bonds, removes this degeneracy. This can be understood 
from the fact that such conformational disorder reduces the overlap between the electron 
orbitals of neighboring carbon atoms and thus increases the chain energy. The energy 
increase is less for the long bonds (with relatively small hopping amplitudes) than for the 
stronger short bonds. 

Due to the random nature of the disorder fluctuations, the preferable sign of the 
dimerization varies along the chain. This explains why domain walls (kinks), separat- 
ing regions with positive and negative dimerization, can be stabilized by disorder .EH 
While in the absence of disorder, kinks (or solitons) are topological excitations with a 
rather high energy (~ 0.5eV in the case of trans-polyacetylene), they do appear in the 
minimal-energy lattice configuration of a disordered chain. This was first noted in nu- 
merical simulations of the Su-Schrieffer-Heerger (SSH) model.10 The large kink energy is 
compensated by adjusting the sign of the chain dimerization to the disorder fluctuations 
in the intervals between the kinks. The weaker is the disorder, the longer the distance 
between neighboring domain walls has to be. Simple arguments, based on an estimate of 
the size of the typical disorder fluctuation that stabilizes a kink-antikink pair, show that 
at weak disorder the average density of disorder-induced kinks in the minimal-energy 
lattice configuration is proportional to the strength of the disorder .0 

The effect of off-diagonal disorder on a Peierls chain is similar to the effect of nonzero 
temperature. It is well-known, that the Z2 symmetry in one- dimensional systems can 
be spontaneously broken only at T = 0. At any non-zero temperature the symmetry is 
restored by thermally induced kinks. In this case the kink creation energy is compensated 
by its large entropy, as kinks can be located at any place in the chain. The topological 
nature of kinks is responsible for the destruction of the long range order in isolated Peierls 
chains, both at arbitrarily small temperature and at arbitrarily weak disorder. 

In this paper we put a firm basis under the results which we obtained through simple 
scaling arguments in Ref.[ |J, by giving a detailed calculation of the free energy and the 
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density of kinks for half-filled Peierls chains with off-diagonal disorder. The continuum 
model describing such chains is introduced in Sec. [TT]. In Sec. f[TT] we briefly repeat the 
arguments of Ref . [ |I| , showing that the energy of a disordered chain can be decreased by 
creating a kink-antikink pair. This implies a special role of kinks. In Sec. [TV] we consider 
the partition function of a disordered chain, treating the lattice classically. Integrating 
out small lattice fluctuations, we obtain an effective free energy describing kinks. The 
details of this integration can be found in the Appendix. We then use the transfer matrix 
approach to reduce the averaging of the chain's free energy over the disorder realizations 
to the averaging of the wave function that describes the relaxation of a spin 1/2 in a 
magnetic field with one random component. For long chains, the latter average can be 
calculated analytically. The derivative of the thus obtained free energy with respect to 
the chemical potential of kinks gives the average density of kinks, induced both thermally 
and by disorder (Sec. |V|). 

As we shall see, in the continuum model with white Gaussian disorder and a classical 
lattice, the entropy of kinks becomes negative below a certain temperature To which 
depends on the disorder strength. Therefore, in Sec. [VI| we also study the generation of 
kinks by disorder in a discrete model, which does not suffer from this pathology. The 
discrete model is the one-dimensional Random Field Ising Model (RFIM). It was realized 
long ago by Imry and Ma,! that kinks destroy the long range order in this model even 
at zero temperature. We present results of numerical simulations of the average density 
of kinks in the RFIM as a function of both the disorder strength and temperature and 



compare to the analytical results for the continuum model. In Sec. [VI I| we summarize and 
conclude. We also connect to previous work and discuss the effects of disorder-induced 
kinks on the optical and magnetic properties of quasi-one-dimensional Peierls systems. 
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II. CONTINUUM MODEL OF A DISORDERED PEIERLS CHAIN 



We start by considering a tight binding model that describes the hopping of electrons 
along a chain of atoms. The electron hopping amplitudes depend on the interatomic 
distances and the relative orientation of the electronic orbitals on neighboring atoms. 
Therefore, the hopping amplitudes are affected by both the lattice motion (the displace- 
ment of the atoms parallel to the chain) and conformational disorder (chain twists). Let 
to denote the hopping amplitude between neighboring atoms in a perfect rigid chain of 
equidistant atoms with lattice constant a. Then, in the presence of atomic displacements 
and conformational disorder, the hopping amplitudes may be written, 

t m ,m+l = to + Cn(u m — U m+ i) + <5t m!m+ i . (1) 

Here, the second term is the SSH-type of electron-phonon interaction,! with the cou- 
pling constant a and u m the displacement of the m'th atom from its uniform-lattice 
position. The third term is a random contribution resulting from the conformational 
disorder. While the lattice displacements u m are dynamic variables, we will assume that 
the fluctuations 8t m , m +i are frozen ("quenched" disorder). 

The Peierls order parameter is the alternating part of the hopping amplitudes, 

A(2ma) = t% m -\^m ~ ^m,2m+i , (2) 

which consists of two parts, 

A(2ma) = A ia <(2ma) + r)(2ma) . (3) 

The first part is the lattice dimerization, 

Ai at (2ma) = a(u 2m -i - 2u 2m + u 2m +i) , (4) 

which describes the alternating part of the hopping amplitude determined by the shifts 
u m of the atoms and is the usual order parameter of the SSH model. The second term 
in Eq.(Q) describes the disorder, 
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7]{2md) — 8t 2 m-l,2m — <^2m,2m+l • (5) 

We assume the random variations of the hopping amplitudes 5t on different bonds to be 
independent. 

For weak electron-phonon coupling and small disorder we can, in analogy to Ref. [ ||, 
use a continuum description of both electrons and lattice with the order parameter 

A(x) = A tat (z) + 7i(x) . (6) 

Here, Ai at (x) is the continuum analog of Eq.(^) and t](x) is the white noise disorder: 

(r)(x)r)(y)) = A8(x - y) . (7) 

Note, that while the random chain twists always decrease the hopping amplitudes 
(<^m,m+i < 0), r](x) can be both positive and negative, as it is the alternating part 
of the fluctuations. 

The Hamiltonian of the continuum model has the form: 

H[A lat {x), n(x)] = E lat [A lat (x)) + H el [A{x)} . (8) 

The first term is the harmonic lattice energy, 

E lat [A lat (x)} = ^— fdxA lat (x) 2 , (9) 



7t\Vf 

where A = is the dimensionless electron-phonon coupling constant (K is the spring 
constant), vf = 2at is the bare value of the Fermi velocity, and we set h — 1 [cf. Ref.[ ||]. 
In this paper we treat the lattice classically, i.e., we disregard the lattice kinetic energy, 
which is reasonable for chains of sufficiently heavy atoms. It should be noted that for 
irans-polyacetylene, which consists of relatively light CH-groups, quantum lattice effects 
may be rather important .EHl£l 

The electrons in the continuum model are described by, 



ip a (x) 



^ V>2<r fa) 
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where the two amplitudes ipi^x) and ip2a{ x ) correspond to particles moving, respectively, 
to the right and to the left with the (bare) Fermi velocity Vp, and a is the spin projection. 
The electron Hamiltonian has the form, 



(10) 



The first term describes the motion of electrons in the presence of both the chain distor- 
tion and the disorder {a% and 03 are the Pauli matrices), while the second term describes 
the (Coulomb) interactions between electrons. Apart from the disorder and the electron- 
electron interaction term the Hamiltonian of our model is the same as the Hamiltonian 
of the continuum version of the SSH model.i 

In the absence of disorder the (half-filled) chain reaches its minimal energy in either 
one of two uniformly dimerized configurations Ai at (x) = ±Ao- This Peierls instability was 
first found for noninteracting electrons^ (H e i_ e i = 0). It also occurs, however, in the pres- 
ence of electron-electron interactions, like, e.g., the on-site Hubbard repulsion U, which 
opens a gap for charge excitations, but (in the absence of electron-phonon interaction) 
leaves the spectrum of spin excitations gapless.0^0 Numerical calculations have shown 
that a moderately large U can even increase the value of the lattice dimerization.&0 As 
we shall shortly see, the existence of two degenerate ground states, which also implies the 
existence of kink solutions, is crucial for the appearance of disorder-induced kinks. The 
precise form of i? e «-eZ does not affect this basic phenomenon; it only is important to the 
extent that it determines the value of Ao and the kink energy. Therefore, in this paper 
we do not specify the expression for H e i-ei- 

An important property of the dimerized Peierls state is the existence of gaps in the 
spectra of spin and charge excitations. For free electrons {H e i- e \ = 0) both gaps are 
equal, while in the presence of Coulomb repulsion the spin gap is smaller than the charge 
gapJli0 In what follows, we will assume the temperature to be much smaller than the 
spin gap, so that we can neglect electronic excitations and replace H e i[A(x)] by its ground 
state expectation value (adiabatic approximation), 
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E el [A{x)} = (0\H el [A(x)]\0) 



(11) 



III. STABILIZATION OF KINKS BY OFF-DIAGONAL DISORDER 

For the sake of completeness, we briefly repeat in this section the arguments of Ref . [ 
|TJ, which show that the lowest-energy lattice configuration in the presence of disorder 
may contain kinks. 

At zero temperature the lattice configuration, i.e., Ai at (x), has to be found by mini- 
mizing the total chain energy, 

E[A lat (x),r](x)} = E el [A{x)\ + E lat [A lat (x)} , (12) 

with respect to Ai at (x) at a given disorder realization r](x). This makes Ai at (x) implicitly 
dependent on r](x). 

As noted in Sec. [Q|, in the absence of disorder the total energy of a half-filled chain 
has two minima, A(x) = ±A , corresponding to two uniformly dimerized configurations 
with the same energy. Apart from the minima, there are infinitely many lattice config- 
urations that are nearly perfect extrema of the total energy. These are the multikink 
configurations, in which a sequence of solitons and antisolitons interpolate between — Ao 
and +Ao and vice versa (Fig. 1). A kink is locally stable, i.e., the chain energy increases, 
when its form is perturbed. The energy of a multikink configuration can be decreased 
only by changing the distances between the kinks. However, when the separation between 
neighboring kinks is large compared to their size (which is of the order of the correlation 
length £ — up/ Aq)^ the change of the energy caused by shifts of the kinks is exponen- 
tially small, so that the energy of the configuration with N kinks is approximately, 

E N = E + NpL . (13) 

Here, [i is the energy needed to create a single kink (about 0.5eV^ for trans-polyacetylene) 
and Eq is the chain energy in the absence of kinks. The kinks can be either charged and 
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spinless, or neutral with spin |.i In the SSH model (H e i- e i = 0), both types have the same 
energy;! if the Coulomb repulsion between electrons is taken into account, the neutral 
soliton has lower energySJii Since the on-site Coulomb repulsion in trans-polyacetylene 
is appreciable,^ we will assume that only neutral kinks are induced by disorder. 

Next we find the change of the energy of a multikink configuration due to weak 
disorder. We will denote the lattice configuration containing N kinks, whose positions 
are described by the N- dimensional vector z = (z±, z 2 , ■ ■ ■ , Zjv)> by An(x\z). To first 
order in r)(x), the correction to the energy of the configuration reads, 

2 r 

5E N = — dxA N (x\z)r](x) , (14) 



ttXvf 

where the extremum condition for the configuration An(x\z) at zero disorder was used. 

For instance, the change of the energy of the uniformly dimerized configuration 
(A(x) = A ) due to disorder is 

2 A r 

5E = — dxr](x) , (15) 



ttXvf 

while for the configuration with an antikink at z\ and a kink at Z2, such that the whole 
disorder fluctuation lies between z\ and Z2, the change of energy equals 5E 2 = ~SE , 
because between z\ and z 2 A 2 {x\z\, z 2 ) ~ — A . We thus see that in the disordered chain, 
the configuration obtained by the perturbation of a kink-antikink pair is energetically 
favorable to the perturbed uniform configuration if 

dxr](x) > 'jkXvp ■ (16) 

Here we introduced 7^ = 7r/x/(2A ), which for free electrons (SSH model) equals 1. Since 
the fluctuations in ^dxri{x) grow with the chain size L, for a sufficiently long chain the 
inequality Eq. (|i~6D can certainly be fulfilled, no matter how small the disorder is. 

We note that in addition to multikink configurations, also non-topological lattice con- 
figurations exist that are extrema of the total chain energy. An example is the polaron.0'§l 
However, since the polaron disturbs the lattice only locally, in an interval of length I ~ £ 
(where £ is the correlation length), a disorder fluctuation having large magnitude in this 



interval is needed to stabilize it. At weak disorder the probability density of such a fluc- 
tuation is exponentially small. On the other hand, the disorder fluctuations stabilizing 
kinks, i.e., those that satisfy Eq. flT6|) , can have small magnitude compensated by a large 
spatial extension given by the distance between neighboring kinks. The weight of such 
a fluctuation grows with its length and becomes of the order unity if the length equals 
the average distance between the disorder-induced kinks.0 In other words, the fluctua- 
tions that stabilize kinks are not at all suppressed. This is why in what follows we take 
into account only multikink configurations, disregarding all other large variations of the 
Peierls order parameter. 

IV. FREE ENERGY OF DISORDERED CHAINS 

We start by considering the partition function of a weakly disordered Peierls chain. 
As was mentioned above, we treat the lattice classically and assume the temperature to 
be sufficiently small to neglect electron excitations across the gap. The partition function 
is then given by the weighted sum over all possible lattice configurations, 

Z[rj(x)} = J DA lat ex V {-(3E[A lat (x), V (x)}} , (17) 

and is a functional of the disorder realization r)(x). In Eq . (|T7|) we have omitted the part 
related to the kinetic energy of the lattice, as ultimately we will only be interested in the 
density of kinks in the chain. 

It is, of course, impossible to perform the integration over Ai at (x) in the partition 
function exactly, since even the dependence of the electron energy on A(x) is, in general, 
not known. However, when the temperature tends to zero, only the lattice configurations 
near the absolute minimum of the chain energy contribute significantly to the partition 
function. As we discussed in the previous section, in the presence of weak disorder the 
minimal energy lattice configuration is close to a certain multikink configuration. How 
many kinks it contains and where they are located, depends on r](x). To make sure that 
for any disorder realization we take all important lattice configurations into account, we 
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perform the integration in Eq.fllTD over A(x) close to all possible multikink configurations. 
Thus, we first perform the saddle-point integration near the configuration An{x\z) with 
fixed number and positions of kinks, and then we integrate over the positions Z\ < Z2 < 
. . . < zn of the kinks and sum over their number N. 

The details of the saddle-point calculation are contained in the Appendix. The re- 
sulting expression for the free energy of a chain with disorder realization 7] (x) is: 



F[r)(x)} = F + F 1 [r 1 {x)} 
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where the first term is the free energy of a chain without kinks and at zero disorder, and 
the second term is given by 

1 



In 



(19) 



Here d(T) is a variable having dimension of length defined by Eq.( A17 ). 

We assume the size of the kinks to be much smaller than the average distance between 
them. In this case the kinks can be replaced by abrupt steps, 

N 



Ajv(a;|z) -> A f N (x\z) = A ]J sign(z n - x) . 



(20) 



n=l 



Equation (p!9|) can then be rewritten as 



F 1 [Tj(x)] = --hi 



where we introduced 



E^/^expj-^^-,/^^^ 



(21) 



(22) 



2A 

7(~\vf 

In the absence of disorder (r](x) = 0) the integration over the kink positions in Eq.(f2~T|) 
is trivial and the free energy can be easily found, 

L 



F l [ V (x) = 0] 



d(T)/3 



(23) 



In this case there are only thermally-induced kinks and their density, 
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is exponentially small at T = f3~ l <C p. This condition, which was already assumed to 
hold in writing Eq . (PD , is always met for solitons in trans-polyacetylene (p ~ 0.5eV). 

Disorder breaks the translational invariance. Nevertheless, the integration over the 
positions of kinks and the sum over their number, can be performed immediately, by 



noting that Eq.(21) can be written in terms of the matrix elements of an ordered expo- 
nential, 



F i[v{ x )] = -~a 1x1 YW\ Tex P 

P „ a' 



f dx (C(^)0"3 + pCTl] 

Jo 



v) , (25) 



where ((x) = /3ki](x), p = e ^ M /rf(T), and \v) is an eigenvector of the Pauli matrix 03, 

a 3 \u) = v\v) . 

The sum in Eq.(^J) over u, v' = ±1 corresponds to the sum over all possible signs of 
the dimerization at the chain ends. We stress that this choice of boundary conditions is 
a matter of convenience and does not affect our final result for the density of solitons. 
The ordered exponential describes the evolution (in imaginary time x) of a spin 1/2 
in a magnetic field, which has a constant x-component and a random time-dependent 
z-component. The corresponding time- dependent Hamiltonian is, 

H(x) = ((x)a 3 + pa 1 . (26) 

Kinks induced in a disordered chain correspond to spin flips and the sum over all possible 
numbers of kinks and the integration over their positions in Eq.(|2~I|) correspond to the 
expansion of the evolution operator in Eq.([2"5D in powers of the random magnetic field. 

It is now convenient to rotate the "coordinate" frame by an angle 7r/2 around e 2 , 
which transforms the Hamiltonian H into 

7T 7T 

H'(x) = exp(i-a 2 )H exp(-i-a 2 ) = -((x)a 1 + pa 3 . (27) 

The expression for the free energy then becomes, 
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F^rjix)] = --In (2(+|Texp 



[ L dxH'(x) 
Jo 



+>J=-iln(2^(L)) 



(28) 



where we introduced the wave function of the spin, ip(x) 
equation, 



, that satisfies the 



with the initial condition ^(0) 

We now introduce v(x) and <b(x), such that 



\ / 



ip = v 



f cosh (f ) 



V sinh (I) J 



The equation for t>(x) 



-— = pcosh(0) 
v ax 



can be easily integrated, yielding 



ipl(L) = cosh 



>(L) 



J exp 


"/ 




Jo 



p / dx cosh(0(x)) 



The equation for 0(x) has the form of a Langevin equation^ 

# dZ7(0) 



dx 



+ /(*) , 



with the potential, 



C/(0) = 2pcosh(0) , 



and the random force, 



f(x) = -2({x) = 2f3K V (x) . 



(29) 



(30) 



(31) 



(32) 



(33) 



(34) 



(35) 



Since the potential U(<fi) grows at — > ±oo, the distribution of at large x tends to 
the equilibrium distribution, given by the Boltzmann formula, 
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P e9 (0) =e T f / d<j>e T f . (36) 

/ J — oo 

where the effective temperature, Tf = 2A(j3n) 2 , is determined by the correlation function 
of the random force, 

(f(x)f(y)} = AA((3k) 2 6(x -y) = 2T f 5(x - y) . (37) 

Assuming the chain length L to be much larger than the average distance between 
kinks, we can use the equilibrium distribution function Eq. (|3~6p to average the free energy 
over the disorder. Using Eqs. (|28|) and (p2|) and neglecting 1/L terms, we obtain, 



(F) = F o -^(cosh(0)) eg . (3* 



The average of cosh(0) over the equilibrium distribution function Eq.(|36|) can be ex- 
pressed in terms of the modified Bessel functions, 



+oo 



KM 



(cosh(0)) eg = / rf0P eg (0)cosh(0) = (39) 



oo 



with z = 2p/Tf oc exp (— (3p) . As we assumed before already that T <^ //, we can use the 
approximate expressions for the modified Bessel functions at small value of the argument 
z. The average free energy of a disordered Peierls chain can then be written in the form, 

(F) =F -LA( -^V 1 . , , (40) 

^ ' \7rXvpJ yU+ |Tln(fi) 1 ^ 

where the temperature T is defined by, 



2 



eXvF \7r7vcr 

and 7 = 1.781072..., is the exponential of Euler's constant. 



1 / AAA 2 , Y 
T o = ~7- — -7= . ( 41 ) 



V. AVERAGE DENSITY OF KINKS 



We now turn to the average density n of disorder-induced kinks, which is obtained by 
differentiating the free energy Eq.fliOD with respect to /x (cf. Eq. (p4]) ). This leads to: 
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which at zero temperature reduces to 

n(T = 0) = - f 2 . (43) 

This zero-temperature result was derived in Ref. [ |IJ] by using scaling arguments. 

Strictly speaking, however, we cannot put T = in our general result, because for 
T < T Q the entropy of kinks, 

(S kink ) = oc In , (44) 



dT VT , 

becomes negative and the density of kinks increases with decreasing temperature. We 
note that in the expression for the free energy Eq.(^), the kink energy /z is effectively 
replaced by, 

H(T) = fx + j-Tln (-j^-) ■ (45) 



2 V T 

This can be readily understood as follows. Each kink enters in the expression Eq. 
with a weight: 

e -MO - l SH e -^ (AQ) 
" d(T) ' 1 } 

where l(T) is the typical size of the thermal fluctuation of the kink position. For the 

white noise disorder that we consider, l(T) can be estimated from the condition, 

{(3n) 2 Al(T) ~ 1 . (47) 

This derives from the criterion that shifting the kink over /(T) from its optimal position, 
should cause a fluctuation of the order of T in its interaction energy Eq.(|14"D with the 
disorder. Using Eq.( |A17|) we find that, 

fH< (48) 

which in combination with Eq.fl46|) indeed gives Eq.(|45|) for yu(T). The temperature 
dependence of /i(T) is plotted in Fig. 2. Below T = T , > 0, i.e., kinks become 
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less "heavy" as temperature decreases, which explains the pathological behavior of the 
entropy. 

Keeping in mind that our description of the lattice is both classical and continuum, 
this pathology is hardly suprising. In practice, quantum effects will prevent the entropy 
from becoming negative at low temperatures. A simpler way to regularize the model, 
lies in a discretization. From Eqs.(f45l) and (46), it follows that To is determined by the 



condition l(T) ~ d(T). For smaller temperatures, l(T) gets smaller than d(T), which 
plays the role of an effective lattice constant and becomes arbitrarily small at small 
temperatures. Clearly, however, the continuum description of the off-diagonal disorder 
fails when 1{T) drops below the lattice constant a. Moreover, the Peierls-Nabarro barrierEEl 
at low temperature fixes the kink position within the unit cell of the atomic lattice. 
These facts motivate us to consider in the next section a discrete version of Eq. (p~8|) , in 
which a temperature-independent lattice constant d is used and the integrations over 
kink positions are substituted by finite lattice sums. 

It should be noted that just replacing d(T) in Eq.flTBD by a temperature-independent 
distance d does not, by itself, make the entropy positive. Retracing the steps following 
Eq.(|T8D in that case leads to 

/ 2A \ 2 1 

(F) = F -LA(— -M ™, (49) 

\ir\v F J /i + 2 Tln(^ 1 ) 



with 



7re\vp \ 7 / 

The density of disorder-induced kinks is then found to be 



A 1 

n — (51) 

(7 fc A^) 2 (i + f ln(^)) 2 ' 



Obviously, this result has the same zero-temperature limit as Eq. (|42|) and it also 
shows a similar pathology below a certain temperature, which now is given by Tq. Only 
if we discretize the kink positions on the lattice with constant d, we will indeed obtain a 
classical model that also at low temperatures gives sensible results. 
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VI. DISCRETE MODEL AND NUMERICAL RESULTS 



As motivated in the previous section, we now consider a discrete model, in which the 
chain is divided into M cells of length d (independent of T) and kinks are allowed to 
be located only at the boundaries x m = md (m = 1, . . . , M) between the cells. Then 
the function f(x) = [cf. Eq.(^)] has a constant value ±1 inside each cell. In other 
words, the chain configuration is described by a set of Ising variables, {cr 1 , cr 2 , . . . , ctm}, 
where a m is the value of f(x) in the mth cell. We impose periodic boundary conditions 
<Jm+i = Kinks occur between cells with opposite values of a and the sum over all 
possible configurations of the Ising spins represents the sum over the number of kinks, 
as well as over their positions on the one- dimensional lattice. This lattice model for the 
statistics of disorder-induced kinks is equivalent to the one-dimensional Random Field 
Ising Modeli (RFIM). The constant-d continuum model which led to Eq.(|5T|) is, in fact, 
the continuum version of this RFIM. 

The partition function of the one-dimensional RFIM reads 

Z = e~P E{um} , (52) 

{tr m } 

with 

M 

E{a m } = Yj 

m=l 

The first term in Eq. (|53"D equals [i times the number of spin flips in the chain and thus 
represents the total kink creation energy. The second term describes the interaction with 
the random "magnetic" field defined by, 

h m = k I dxr/(x) , (54) 

i.e., the value of the field in each cell is proportional to the value of disorder in the 
continuum model integrated over the cell. The distribution of the random field is also 
Gaussian, with correlator: 

(h k h t ) = DS kl , (55) 
16 



[\ — a m a 



m ^ m+1 , 



(53) 



where D = dn 2 A. 

It is well-known that arbitrarily small disorder or temperature destroy the long range 
order in the one-dimensional RFIM.i The density of the thermally excited domain walls 
in the absence of disorder is e - ^, while at zero temperature the density of the disorder- 
induced kinks is proportional to the strength of the disorder. A full calculation of the 
average free energy and the density of kinks in the RFIM involves numerical simulations. 
We performed such simulations using the transfer-matrix approach, as we did in the con- 
tinuum treatment. In the discrete case, the Schrodinger equation (|29|) for the continuum 
model is replaced by 



VWi = T m ip m , (56) 



with the transfer matrix T m , 



1 exp{[3h m } exp{-/5(/i - h m )} ^ 



T 



(57) 



y exp{-/5(/i + h m )} exp{-(3h m } j 
The free energy for a given realization {h m } of the random field is then given by 



F[{h m }] = -T In (tt f[ f m ) . 

\ m=l J 



(5f 



To obtain a smooth temperature dependence of the density of kinks, we had to average 
the free energy over 10 3 random field realizations for a chain with 10 3 sites. 

We now turn to the results of these simulations. In Fig. 3 the free energy is plotted as 
a function of temperature for two values of the disorder strength: D = 0.005/i 2 and D = 
O.OlyU 2 . Stars indicate the numerical results for the discrete model, while the solid lines 
represent the continuum version of the RFIM (Eq. (|49D ). The latter curves have maxima 
at T = Tq, below which the entropy becomes negative. For convenience, we also plotted 
the dashed curve, which passes through the maxima of the solid lines. The continuum 
results become meaningless to the left of this curve. The slope of the numerically obtained 
free energy of the discrete model is always negative, corresponding to a positive entropy, 
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which tends to as T — ► 0. We also observe that above Tq the numerical (discrete) and 
analytical (continuum) results only differ little from each other. 

In Fig. 4 we give the density of kinks as a function of temperature for the same two 
values of the disorder as in Fig. 3. Again, solid lines and stars represent the continuum 
model and the discrete model, respectively, while the dashed curve passes through the 
minima of the continuum results, below which the continuum model becomes meaningless. 
The temperature dependence of the density of kinks for T < 0.15// is very slow. At 
T = 0.15// the density of thermally induced kinks in a chain without disorder is about 
1.3 • 10~ 3 /d. This means that the increase of the kink density observed in Fig. 4 is mostly 
related to the fact that more multikink configurations become available in a disordered 
chain as temperature grows. 

Figure 5 shows the density of kinks as a function of the disorder strength D at T = 0. 
Stars give the numerical results for the discrete model (simulated at T = 0.005// in order 
to avoid numerical instabilities). The dashed line gives the T = continuum result, 
Eq.(f4*3"D, which in terms of D reads n(T = 0) = D/(dfi 2 ). As is observed from Fig. 4, for 
decreasing disorder these two results get in better agreement with each other. This is a 
consequence of the fact that Tq then also decreases. In the limit of D — > 0, the agreement 
should become exact, as indeed appears from Fig. 5. For higher disorder, however, 
the nonphysical nature of the continuum model for T < Tq causes the numerical and 
analytical results to deviate by an increasing amount at T = 0. An alternative approach 
to obtain an analytical zero temperature result, is to use the continuum result Eq . fl5l"l) 
at the lowest possible value where the result is still physical, i.e., at T = Tq. The thus 
obtained value of the density of kinks as a function of the disorder is also plotted in Fig. 5 
(solid line) and is indeed seen to be in much better agreement with the numerical T = 
results over a large interval of the disorder. This shows the importance of the corrections 
to the linear D dependence, expected on the basis of scaling arguments,!! of the density 
of domain walls in the RFIM. 
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VII. SUMMARY AND DISCUSSION 



In this paper we studied a continuum model describing Peierls chains with a doubly 
degenerate ground state in the presence of weak off-diagonal disorder. In such chains, 
thermally- and disorder-induced neutral solitons (dimerization kinks) occur. By inte- 
grating out all irrelevant lattice degrees of freedom, we obtained the free energy of these 
kinks. This free energy was then averaged over the disorder realizations leading to an 
analytical expression Eq.(ff2|) for the average density of kinks. 

In the continuum model, the entropy of the kinks was found to become negative 
below a certain temperature T . At this temperature the thermal fluctuation of the 
kink positions l(T) becomes of the order of the effective "lattice constant" d(T). It was 
then argued that below this temperature one should either take the effects of the lattice 
discreteness into account, or consider the lattice quantum-mechanically. We noted that 
the discrete version of the model that describes the thermal or disorder-induced kinks, is 
the one-dimensional RFIM. The free energy and the density of kinks in the latter model 
were found numerically and the results were compared with the continuum model. It was 
found that for T > T , both models are in good agreement. 

To obtain the free energy of the kinks (Appendix A), we used a technique developed 
for the semiclassical solution of quantum problems.il This is not accidental, because a 
Peierls chain with a doubly degenerate ground state in the presence of weak disorder 
or at a small temperature behaves similarily to a quantum particle in a double well. If 
the wells are very deep, the small quantum fluctuations around the two classical vacua 
may be relatively insignificant. A more important effect, however, is that the quantum 
particle can tunnel between the two minima. The tunneling probability is described by 
instantons, which are the imaginary-time classical trajectories leading from one minimum 
of the potential to the other. In fact, an exact mapping exists between the instanton 
gas for a double-well potential and the gas of thermally-induced kinks.il The disorder- 
induced kinks are somewhat different: their positions, unlike the instanton positions, are 
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not arbitrary. These kinks are, to some extent, pinned by the disorder which induces 
them. 

The main assumption of our calculation is the validity of the expansion of the free 
energy in powers of the disorder. As was already discussed in Ref. [ [TJ, this translates into 
the condition: 

<o « 1 , (59) 

meaning that the average distance between kinks should be much larger than their size. 
For stronger disorder the lattice distortions would be random over the length scale of the 
order of £o, and the notion of multikink configurations would lose sense. 

Effects of disorder on the properties of Peierls chains have been studied before by a 
number of authors. Many of these studies were devoted to charged solitons induced by 
random doping! and the effects of localized bond and site impurities on the electronic 
states.il The minimal-energy lattice configuration in the presence of a finite density 
of bond impurities was studied numerically in Refs. [ |!|f§]. There it was found that 
for some disorder realizations the configuration contained kinks. As these calculations 
were performed within the SSH model (in which electron-electron interactions are not 
included), the obtained disorder- induced kinks did not have a definite charge or spin. 
Also, no attempt was made in these papers to obtain (numerically or analytically) the 
density of solitons as a function of the disorder strength. 

Another large body of work on disordered Peierls systems deals with the Fluctuating 
Gap Model (FGM).S This work not only addresses the study of disorder per se,00 but 
also includes studies of quantum lattice fluctuations modeled by static disorder.&B In 
the FGM, one assumes that the lattice part A/ at of the order parameter is a constant 
throughout the chain and does not respond to the local disorder. This immediately elimi- 
nates the possibility of disorder- induced solitons. Nevertheless, solitons play an important 
role in the FGM: we recently showed that the optimal fluctuation that determines the 
density of states inside the pseudogap in this model has the form of a soliton-antisoliton 



pairH 
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We finally mention that, while we considered weak off-diagonal disorder, one may 
also consider a type of disorder, where the electron conjugation over certain bonds is 
broken completely (£ n ,n+i = 0). This model leads to a distribution of chain sizes.S If 
the conjugation breaks occur randomly, chains with an odd number of carbon atoms will 
occur, which necessarily contain a kink.0 Thus, the density of kinks is proportional to 
the disorder strength (i.e. the density of chain breaks), as it is in our model (see Eq. fl43f) ). 

As was already discussed in our previous paper,0 the electron state in the presence 
of neutral solitons is different from that of a perfectly dimerized chain. This gives the 
possibility to detect the disorder-induced kinks. First, the kinks should give rise to 
a peak in the optical absorption spectrum that occurs inside the gap of the perfectly 
dimerized chain. Second, having spin 1/2, the neutral kinks contribute to the Curie 
susceptibility. As we noted in Ref.[ Q, experiments on undoped trans-polyacetylene do 
not seem to agree with the large density of free solitons that one would expect on the basis 
of the small average electron conjugation length (several tens of carbon atomsS) in this 
material. In particular, experiments indicate that only 1 free spin occurs per 3000 carbon 
atoms.il It should be noted, however, that a meaningful comparison between theory and 
experiment can only be made if one accounts for interchain interactions, which suppress 
the generation of isolated kinks.il A study of the interplay between the disorder-induced 
kinks, destroying the long range coherence, and the interchain interactions, binding kinks 
into pairs, will be reported elsewhere. It should be noted that this interplay seems to be 
very important to describe the Peierls phase transition.0 

We finally note that another important issue is the quantum treatment of the lattice 
dynamics coupled to the spin dynamics of the kinks. The interaction between neighboring 
kinks may bind their spins into a singlet ground state, thus explaining why the density 
of free spins observed in undoped trans-polyacetylene is so low. We plan to address this 
problem in the future. 
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APPENDIX A: 

In this Appendix we obtain the effective free energy of kinks by integrating out all 
irrelevant lattice degrees of freedom in the partition function of the disordered chain 
Eq.(0). 

It is convenient to replace the integration over lattice configurations A/ at by the inte- 
gration over A (a;), which also includes the disorder [cf. Eq.(^)] . The chain energy can 
be written as, 

E[A lat (x),r](x)} = E el [A(x)} + E lat [A(x) - r)(x)} . (Al) 
Since the lattice energy is quadratic in A/ at = A(x) — T](x), we obtain, 

E[A lat (x), rj(x)) = E[A(x), 0] - — ^— f L dxA(x)r](x) + — — fax^x) 2 , (A2) 

TVAVp JO TlAVp Jo 

where the first term is the chain energy at zero disorder and L is the chain length. 

For A(x) close to the multikink configuration An(x\z) we can expand the chain energy 
E[A(x), 0] in powers of £(x) defined by: 

£(x) = A(x) - A N (x\z) . (A3) 

Retaining terms up to second order in £(x), we have 

E[A, 0} = E + Nfi + — — / dx£(x)D N (z)£(x) . (A4) 

n\v F Jo 

Here, we assumed the distances between kinks to be much larger than the correlation 
length £o- As we discussed in Sec. [II], in the absence of disorder the multikink configu- 
ration is then an almost exact extremum of the chain energy, which is why Eq. (|A4|) does 
not contain a first-order term. The operator Dn( z ) is defined through: 
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D N (z)£(x) = JdyD N (x,y\z)£(y) , (A5) 

with the kernel Djs[(x,y\z) related to the second variational derivative of the total chain 
energy: 



1 D (x xlz)- 1 52E[A '° ] 
-JJ N {Xi, x 2 \z) 



(A6) 

A=A N (x\z) 



7t\v f v ' 1 ' 2 5A(x 1 )5A(x 2 ) 

Let Q = Jff denote the bare optical phonon frequency. Then, the eigenf unctions <p a (x\z) 
of D N (z) obeying, 

D N (z)cf> a (x\z) = {^-) *a{x\*) , (A7) 

are (optical) phonon modes in the presence of kinks with frequencies w a (z). 

When the kinks are far apart, the A^ lowest phonon frequencies (a = 1, . . . , N) are 
very small, since these normal modes correspond to shifts of the kinks, which leave the 
chain energy practically unchanged. (The frequency of one mode is exactly zero, as 
it corresponds to the translation of the multikink configuration as a whole). To avoid 
double counting, the integration over the first A^ modes is replaced by the integration over 
the A^ kink positions, which is the standard approach in the instanton calculus.@ Then 
the measure of the functional integration in the vicinity of the multikink configuration 
becomes, 

D£(x) = J N d N z J] d£ a . (A8) 

a>N 

where £ Q are the amplitudes of the normal modes, 

^HE^W 2 )' ( A9 ) 

a 

and Jn is the Jacobian of the substitution of the integration over the A^ (almost) zero 
modes by the integration over the kink positions. We have: 

J N = {J X ) N = cf , (A10) 

with Ji the Jacobian for the single- kink configuration Ai(x), leading to 
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The constant c can be expressed as c = tt\vfQ 2 Mi s /2, with the soliton mass, which 
in the SSH model is about 6m e Jl 

The Gaussian integration over the amplitudes of the non-zero modes now reduces the 
partition function to, 



^ = ^oE(^!r) 2 fd N zR N (z)exp{-[3W N (z)} . (A12) 



N 



N x n 2 Xv F 

Here 

a>N w a\ 

with r] a the expansion coefficients of rj(x) over the orthonormal basis <f) a (x\z) (cf. 
Eq.([MS|)). Z is the partition function of the chain in the absence of kinks and disor- 



W N (z) = Nfi- [ L dxA N (x\z) v (x) + — !— ( [ L dx V (x) 2 - -rr^vl ) , (A13) 

nXv F Jo n\v F \ Jo XirUiliz) I 



der: 



z --1Wts) ■ (A14) 

where w a (0) denote the phonon frequencies in the absence of kinks and disorder. Finally, 
the factor -Rat(z) in Eq.( |A12| ) is the product of nonzero phonon frequencies for the lattice 
configuration containing no kinks divided by the same product for the configuration with 
N kinks, 

^)^(n^)Cn^)" 1 - 

Similarily to Eq. ( |A10| ), for widely separated kinks, i?iv(z) equals the iVth power of the 



regularized product of nonzero phonon frequencies for a configuration with one kink,@ 



R N (z) = (R 1 (z 1 )f = r^ } (A16) 

where r is a number of the order unity. 

Finally, for weak disorder we can neglect the third term in Eq. (|A13|) , as it is of second 



order in rj(x). From Eqs.( |A12| ), (|A13| ), and ( [A16Q we then obtain the expressions Eqs. (|T8|) 
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and fll^D for the free energy of a chain with disorder realization r)(x), where d(T) is given 
by 



«T>=l\l*£~H%- (A17) 

In this calculation we neglected electron excitations, assuming the temperature to be 
sufficiently low. However, the electron ground state for the lattice configuration contain- 
ing N neutral kinks, separated by distances much larger than the correlation length, is 
2 JV -fold degenerate, due to the fact that the spin projection of each kink can be arbitrary. 
Because of this degeneracy, the weight factor of each kink acquires an extra factor of 2, 
which was taken into account in Eq.( A12|) . The degeneracy is not exact: excitations of 
this spin system have a typical energy of the order of the exchange constant for the spins 
of two neighboring kinks. Thus, our assumption requires that the temperature is much 
larger than this energy [A exp(— l/(n£o))], yet much smaller than the Peierls gap [A ]. 
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FIGURE CAPTIONS 



FIG.l Order parameter along a Peierls chain for a multikink configuration with 7 
kinks. Zi, Z2, ■ ■ ■ , Z7, denote the kink positions. 

FIG. 2 Temperature dependence of the effective chemical potential of kinks according 
to Eq.(H) for /1 = 10T . 

FIG. 3. The free energy of dimerization kinks in a Peierls chain as a function of 
temperature for two values of the disorder strength: D = 0.005/i 2 (a) and D = 0.01/z 2 (b). 
The solid lines give the analytical result Eq. (|49|) , while the stars indicate the numerical 
results. The part of the free energy related to an ordered chain without kinks (F ) is 
discarded. M = L/d denotes the number of unit cells. The dashed curve intersects the 
solid lines at the temperature for which the analytically calculated entropy equals zero. 

FIG. 4. The density of dimerization kinks in a Peierls chain as a function of tempera- 
ture for two values of the disorder strength: D = 0.005/i 2 (a) and D = O.Ol/i 2 (b). The 
solid lines give the analytical result Eq. ([51]) and the stars indicate the numerical results. 
Dashed curve as in Fig. 3. 

FIG. 5. Density of dimerization kinks in a Peierls chain as a function of disorder 
strength. The dashed line represents the analytical continuum result Eq. ([)]]) taken at 
T = 0, the solid line gives the same result at T = Tg, while the stars indicate the results 
of numerical simulations of the discrete model at T = 0.005/i (which for all practical 
purposes may be considered zero temperature). 
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FIG. 5 M.V. Mostovoy, etal. 




FIG. 4 M.V. Mostovoy, etal. 
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